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I review recent theoretical advances in quantum chromodynamics. Particular emphasis is put on developments related 
to the precise prediction and interpretation of experimental data from present and future high energy colliders. 



1 Introduction 

Quantum Chromodynamics (QCD) is well estab- 
lished as theory of strong interactions through a large 
number of experimental verifications. The era of 
'testing QCD' is clearly finished, and QCD today 
is becoming precision physics. The next generation 
of high energy collider experiments are all performed 
at hadron colliders, where (in contrast to LEP and 
SLC) QCD is ubiquitous. Any precision measure- 
ment (strong coupling constant, quark masses, elec- 
troweak parameters, parton distributions) at the 
Tevatron and the LHC, as well as any prediction of 
new physics effects and their backgrounds, relies on 
the understanding of QCD effects on the observable 
under consideration. 

The derivation of precise QCD predictions for 
collider observables poses several theoretical and 
computational challenges. The most important chal- 
lenge is the fact that QCD describes quarks and glu- 
ons, while experiments observe hadrons. This mis- 
match is either accounted for by a description of 
the parton to hadron transition through fragmen- 
tation functions or by defining suHiciently inclusive 
final state observables, such as jets. Moreover, the 
strong coupling constant is considerably larger than 
the electromagnetic coupling constant at scales typ- 
ically probed at colhders: as{Mz) — 15acm(-^^z), 
resulting in a slower convergence of the perturbative 
expansion. As a consequence, a precise description of 
QCD observables (precise means here that the theo- 
retical uncertainty becomes similar to the achieved or 
projected experimental errors) is obtained only by in- 
cluding higher order corrections, often requiring be- 
yond the next-to- leading order. The largeness of the 
strong coupling constant also implies that multipar- 
ticle final states are rather frequent. Finally, many 
collider observables involve largely different scales, 
such as quark masses, transverse momenta and vec- 
tor boson masses. These give rise to potentially large 
logarithms, which might spoil the convergence of the 



perturbative series and need to be resummed to all 
orders. 

In this talk, I shall try to highlight recent theo- 
retical progress towards precision QCD at colliders, 
focusing on heavy quark production in Section 2, on 
jets and multiparton final states in Section 3, on pho- 
tons in Section 4 and electroweak bosons in Section 5. 
Finally, a summary of the current state-of-the-art 
and of yet open issues is given in Section 6. 

2 Heavy Quarks 

Heavy quark production is one of the main top- 
ics investigated at high energy collider experiments. 
Heavy quarks are of particular interest to probe the 
fiavour sector of the standard model, which is less 
well tested than the gauge sector. Also, many ap- 
proaches to physics beyond the standard model, of- 
ten related to electroweak symmetry breaking and 
mass generation, predict new effects to be most pro- 
nounced in observables involving heavy quarks. 

2.1 Total Cross Sections 

The total cross sections for the production of heavy 
quarks can be computed reliably within perturba- 
tion theory. The current state-of-the-art is a next- 
to-leading order (NLO) calculation^ , which is further 
improved by summing large logarithms due to soft 
gluon emission up to the next-to-leading (NLL)^ and 
next-to-next-to-leading logarithmic (NNLL)'^ level. 
As can be seen from Figure 1, these predictions are 
in good agreement with experimental data on the 
total tt cross section at the Tevatron* and the total 
bb cross section at HERA-B^ (which both actually 
refer to similar kinematical values of mq/^/s). The 
theoretical uncertainty on the prediction for HERA- 
B is larger for two reasons: the larger value of the 
strong coupling at rub than at mt and the dominance 
of gg initial states in pN collisions (HERA-B) com- 
pared to qq dominance in pp collisions (Tevatron). 
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Figure 1. Total cross sections for tt at the Tevatron^ and bb 
at HERA-B^. 



The effects due to soft gluon resummation turn out 
to be moderate, but do yield a significant decrease in 
the uncertainty of the theoretical prediction. Further 
uncertainties on the prediction of the ti cross section 
come from the parton distribution functions^. 

2.2 Transverse Momentum Distributions 

Differential distributions of hadrons containing b 
quarks measured in hadron-hadron, photon-hadron 
or photon-photon collisions have been in apparent 
discrepancy with theoretical predictions for quite 
some time. The spectrum of hadrons measured 
at CDF'' is one of the most recent examples for this 
discrepancy. 

The theoretical prediction for B meson produc- 
tion involves a convolution of the hard matrix ele- 
ment for heavy quark production in parton-parton 
scattering with initial parton distributions and fi- 
nal state fragmentation functions describing the non- 
perturbative transition from a 6 quark to a _B hadron. 
It is in particular the latter which is suspected to 
account for the discrepancy between theoretical pre- 
diction and experimental data, especially since it has 



been observed^ that the transverse momentum dis- 
tribution of ^-tagged jets^ (which has little sensitivity 
to fragmentation functions) is in much better agree- 
ment with theoretical predictions. 

The definition of heavy quark fragmentation 
functions is not free from ambiguities, since some 
aspects of these functions are actually calculable in 
perturbation theory^''. In extracting these fragmen- 
tation functions from data on B hadron production 
in e+e^ collisions, several choices are made, related 
to the order of perturbation theory, the incorpora- 
tion of mass effects in the matrix elements, the re- 
summation of potentially large perturbative terms, 
the inclusion of power corrections^^, the correction 
of data for parton showers or the parametric form 
of the ansatz used in the determination. Unfortu- 
nately, the sensitivity of the fragmentation function 
on the assumptions used in the extraction from e^e^ 
spectra is often overlooked when using this fragmen- 
tation function to compute heavy hadron spectra at 
colliders. 

Recently, an approach incorporating quark mass 
effects, perturbatively calculable components of the 
heavy quark fragmentation function^*' and resumma- 
tion of large logarithms up to the next-to-leading 
logarithmic level has been put forward with the 
fixed-order next-to- leading log (FONLL) scheme-"^^. 
This approach requires only a small, genuinely non- 
perturbative component of the fragmentation func- 
tion to be fitted to e~^e~ data. In order to expose 
the information content actually relevant to heavy 
hadron spectra at hadron colliders, this fit is done in 
moment space. 

In view of new data from ALEPH^^, a phe- 
nomenological study of B hadron production at col- 
liders based on the FONLL scheme was performed^"*. 
It was shown that the consistent treatment of the 
fragmentation function in extraction and prediction 
reduced the discrepancy between data and theoreti- 
cal prediction considerably. The theoretical predic- 
tion is however still falling somewhat short of the 
experimental data, which is probably due to cur- 
rently uncalculated corrections beyond NLO. More 
recently, the same framework was applied to charmed 
hadron production at hadron collider s^'"^. In this 
case, one also observes that the experimental data^^ 
exceed the theoretical prediction. Figure 2, although 
the effect is less pronounced than for bottom hadron 
production. Comparison of massless^^ and massive^^ 
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Figure 2. Transverse momentum spectrum of _D+ hadrons at 
CDF, compared to calculation using FONLL fragmentation 
functions^^ . 

calculations for hadronic charm production^^ indi- 
cates only moderate mass effects, which are smaller 
than the theoretical uncertainty. 

Many collider experiments also report an excess 
in the b quark production spectra. In interpreting 
these data, it must always be kept in mind that 
it is not b quark but B hadron production which 
is observed in the experiment. Information on b 
quark production is only inferred from these data 
using some model for the heavy quark fragmenta- 
tion. As discussed above, there are numerous am- 
biguities, which can yield inconsistent predictions if 
not implemented consistently. In view of the rather 
sizable effects due to a consistent treatment of the 
fragmentation function observed on the B hadron 
spectra at CDF, it might be that the data sets on 
b quark spectra have to be reanalyzed incorporating 
the new experimental and theoretical information on 
the b quark fragmentation functions in a consistent 
manner. 

3 Jets and Multiparticle Production 

Hadronic jets at large transverse momenta are pro- 
duced very copiously at colliders. Final states with 
a small number of jets are measured to very high 
experimental accuracy, such that they can be used 
for precision measurements of the strong coupling 
constant and of parton distribution functions. For 
these measurements, the uncertainty on the theoret- 
ical prediction is often the dominant source of error. 



and one would consequently like to have theoreti- 
cal calculations to be more accurate, which implies 
in general an extension towards higher perturbative 
orders. Multiparton final states, involving a large 
number of jets, can on the other hand mimic final 
state signatures induced by physics beyond the stan- 
dard model, thus forming an irreducible background 
to searches. For these, QCD predictions serve as 
a guidance to devise search strategies, and one de- 
mands QCD to yield a description of the full hadronic 
final state. 



3.1 Leading order calculations 

Multiparton final states are described using lead- 
ing order QCD predictions, implemented in flexible 
multiparton matrix element generators. These pro- 
grams evaluate the scattering amplitudes using effi- 
cient representations of helicity amplitudes or fully 
numerically from the interaction Lagrangian. Exam- 
ples of these codes are VECBOSl^ COMFHEPi^, 
MADGRAPH2", GRACE^i, HELAC^^, 

ALPHGEN23 and AMEGIC++^^. Using these, the 
computation of 2 — > 8 reactions is feasible on cur- 
rent computers. These programs are then combined 
with automatic integration over multiparticle phase 
spacs, using for example RAMBO^^, PHEGAS^'' or 
MADEVENT^"^. Most programs can be interfaced 
to hadronization models using standard interfaces'^. 

Matrix element calculations accurately include 
large angle single gluon radiation. At small angles 
from the emitting particles, one does however en- 
counter multiple gluon radiation, which can be ac- 
counted for by parton showers. Recently, a generic 
procedure was devised to combine both descriptions 
for multiparton final states in a modified matrix ele- 
ment plus vetoed parton shower'^. 

Including these developments, leading order 
QCD provides the basis of Monte Carlo event gen- 
erators. However, its predictions contain large (and 
non-quantifiable) errors due to the setting of renor- 
malization and factorization scales. Leading order 
QCD is therefore good tool to estimate relative mag- 
nitudes of processes and to design searches. Once 
precision is required (e.g. to identify a discovery with 
a particular model), it is not sufficient. 
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Figure 3. Dependence of the NLO prediction for the 3j cross 
section at the Tevatron on renormalization and factorization 
scale'^-^. 
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Figure 4. Transverse momentum distribution of top quark 
pairs at the LHC as predicted by MCONLC^. 



3.2 Next-to-leading order calculations 

Including ncxt-to-leading order QCD corrections im- 
proves the theoretical predictions in numerous ways 
by reducing the renormalization scale uncertainty, 
providing reliable normalizations of cross sections, 
and reliable error estimates. Moreover, NLO is the 
first order where differences between jet algorithms 
show up. In contrast to leading order, there is no 
generic procedure for doing NLO calculations, such 
that each new process under consideration implies a 
completely new calculation. 

For hadron colliders, NLO results are available 
for all relevant 2^2 reactions; 2^3 reactions 
are the current frontier. A number of 2 — > 3 results 
(each involving several man-years of work) became 
available recently: pp ^ V -\- 2j^°, ep ^ (3 -t- 
pp 3j'^^, pp — *■ pp ttH^'^ and the vector 

boson fusion processes H-{-2j^^, pp V-\-2j^^. 
Some of the features of NLO calculations, such as 
the improved scale dependence and the differences 
between jet algorithms are illustrated in Figure 3. 

To overcome the large amount of work required 
for each NLO calculation, efforts are under way to- 
wards their automatization. The NLO calculation 
for an n parton reactions contains the one-loop n 
parton matrix elements, the tree level n -\- 1 parton 
matrix elements and a procedure to extract the in- 
frared singularities from both and to combine them. 
While this procedure has been automatized for the 



tree level real radiation matrix elements long ago'^ , 
there is at present no automatic procedure to com- 
pute one-loop integrals. Very recently several algo- 
rithms were proposed, including the analytic reduc- 
tion of hexagon integrals'^'*, a subtraction formalism 
for virtual corrections'^^ and the numerical evalua- 
tion of hexagon integrals^°. 

Another important development is the combina- 
tion of NLO calculations with parton showers, as re- 
alized in the MC@NLO approach""^. This approach 
introduces a modified NLO subtraction method, 
where both real and virtual contributions become ini- 
tial conditions for the parton shower. In this, hard 
radiation is accurately described by the NLO matrix 
element, while multiple soft radiation is accounted 
for by the parton shower; a double counting of con- 
tributions is avoided. So far, this formalism has been 
applied to VV, bb and production at hadron colhd- 
ers. Figure 4 illustrates that MC@NLO^^ smoothly 
connects the kinematic region dominated by multiple 
radiation at small transverse momenta to the region 
controlled by single hard radiation at large transverse 
momenta. 

3.3 Next-to-next-to-leading order calculations 

Despite the evidently good agreement of NLO QCD 
with experimental data on jet production rates, pre- 
dictions to this order are insufficient for many appli- 
cations. For example, if one uses data on the single 
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Figure 5. Errors on extraction of as from single jet inclusive 
cross section at CDF . 



jet inclusive cross section'*'^ compared to the NLO 
theoretical prediction'*^'^'^ to determine the strong 
coupling constant as , it turns out that the dominant 
source of error on this extraction is due to unknown 
higher order corrections. Given that the theoretical 
prediction to infinite order in perturbation theory 
should be independent of the choice of renormaliza- 
tion and factorization scale, this error can be esti- 
mated from the variation of the extracted as under 
variation of these scales, as seen in Figure 5. As a 
result, CDF find from their Run I data 

as{Mz) = 0.1178 ±0.0001(stat)t°:^^(sys) 
^0.0047 (scale) ± 0.0059(pdf). 

It can be seen that the statistical error is already neg- 
ligible; improvements in the systematic error can be 
anticipated in the near future. To lower the theoret- 
ical error, it is mandatory to compute next-to-next- 
to-leading order (NNLO) corrections to the single jet 
inclusive cross section. 

A similar picture is true in e~^e~ annihilation 
into three jets and deep inelastic (2 + 1) jet produc- 
tion, where the error on the extraction of as from 
experimentally measured jet shape observables'*^ is 
completely dominated by the theoretical uncertainty 
inherent in the NLO QCD calculations. 

Besides lowering the theoretical error, there is a 
number of other reasons to go beyond NLO in the 
description of jet observables. While jets at NLO 
are modeled theoretically by at most two partons. 



NNLO allows up to three partons in a single jet, 
thus improving the matching of experimental and 
theoretical jet definitions and resolving the internal 
jet structure. At hadron colliders, NNLO does also 
account for double initial state radiation, thus pro- 
viding a perturbative description for the transverse 
momentum of the hard final state. Finally, including 
jet data in a global NNLO fit of parton distribution 
functions, one anticipates a lower error on the pre- 
diction of benchmark processes at colliders. 

The calculation of jet observables at NNLO re- 
quires a number of different ingredients. To com- 
pute the corrections to an n-jet observable, one needs 
the two- loop n parton matrix elements, the one- 
loop n + 1 parton matrix elements and the tree 
level n + 2 parton matrix elements. Since the lat- 
ter two contain infrared singularities due to one or 
two partons becoming theoretically unresolved (soft 
or coUinear), one needs to find one- and two-particle 
subtraction terms, which account for these singu- 
larities in the matrix elements, and are sufficiently 
simple to be integrated analytically over the unre- 
solved phase space. One-particle subtraction at tree 
level is well understood from NLO calculations'^^ and 
general algorithms are available for one-particle sub- 
traction at one loop"*^, in a form that could recently 
be integrated analytically'**. Tree level two-particle 
subtraction terms have been extensively studied in 
the literature'*^, their integration over the unresolved 
phase space was up to now made only in one partic- 
ular infrared subtraction scheme in the calculation 
of higher order corrections to the photon-plus-one- 
jet rate in e~^e~ annihilation^*^. The same techniques 
(and the same scheme) were used very recently in the 
rederivation of the time-like gluon-to-gluon splitting 
function from splitting amplitudes'^. A general two- 
particle subtraction procedure is still lacking at the 
moment, although progress on this is anticipated in 
the near future. 

Concerning virtual two-loop corrections to jet- 
observables related to 2 — > 2 scattering and 1^3 
decay processes, enormous progress has been made 
in the past years. Much of this progress is due to 
several technical developments concerning the eval- 
uation of two-loop multi-leg integrals. Using itera- 
tive algorithms'^^, one can reduce the large number of 
two-loop integrals by means of integration-by-parts^'^ 
and Lorentz invariance^'* identities to a small num- 
ber of master integrals. The master integrals relevant 
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to two-loop jet physics are two-loop four-point func- 
tions with all legs on-shell^^ or one leg off-shell^^, 
which were computed using explicit integration or 
implicitly from their differential equations^^. 

Combing the reduction scheme with the master 
integrals, it is straightforward to compute the two- 
loop matrix elements relevant to jet obscrvablcs us- 
ing computer algebra^^. Following this procedure, 
massless two-loop matrix elements were obtained 
for Bhabha scattering"''^, parton-parton scattering 
into two partons^^, parton-parton scattering into 
two photons®^, as well as light-by-light scattering^^ . 
Two-loop corrections were also computed for the off- 
shell process 7* qqg^"^, relevant to e+e~ 3j. 
Part of these results were already confirmed^'^ using 
an independent method^*. Related to e+e~ 3j by 
analytic continuation^^ are (2-1- l) j production in ep 
collisions and V+j production at hadron colliders. A 
strong check on all these two-loop results is provided 
by the agreement of the singularity structure with 
predictions obtained from an infrared factorization 
formula^^ . 

More recently, first results were obtained for 
master integrals involving massive internal propaga- 
tors, as appearing in the two-loop QED corrections 
to the 7* QQ vertex^'' or in the two-loop elec- 
troweak corrections to the V ^ qq vertex^^. 

4 Photons 

Photons and gauge bosons provide very prominent 
final state signatures at colliders. Their study allows 
the precise determination of electrowcak parameters 
at hadron colliders, and their final state signatures 
are often background to searches, such as photon pair 
production to the Higgs search in the lower mass 
range. 

4-1 Isolated Photons 

Photons produced in hadronic collisions arise essen- 
tially from two different sources: 'direct' or 'prompt' 
photon production via hard partonic processes such 
as qg — > q-y and qq g-y or through the 'fragmen- 
tation' of a hadronic jet into a single photon carry- 
ing a large fraction of the jet energy. The former 
gives rise to pcrturbatively calculable short-distance 
contributions whereas the latter is primarily a long 
distance process which cannot be calculated pertur- 



batively and is described in terms of the quark-to- 
photon fragmentation function. In principle, this 
fragmentation contribution could be suppressed to 
a certain extent by imposing isolation cuts on the 
photon. Commonly used isolation cuts are defined 
by admitting only a maximum amount of hadronic 
energy in a cone of a given radius around the photon. 
An alternative procedure is the democratic clustering 
approach suggested in^^, which applies standard jet 
clustering algorithms to events with final state pho- 
tons, treating the photon like any other hadron in 
the clustering procedure. Isolated photons are then 
defined to be photons carrying more than some large, 
predefined amount of the jet energy. 

Both types of isolation criteria infrared safe, al- 
though the matching of experimental and theoret- 
ical implementations of these criteria is in general 
far from trivial. It was pointed out recently^^ that 
cone-based isolation criteria fail for small cone sizes 
R (once aglnii"^ 1), since the isolated photon 
cross section exceeds the inclusive photon cross sec- 
tion. This problem can only be overcome by a resum- 
mation of the large logarithms induced by the cone 
size parameter. The contribution from photon frag- 
mentation to isolated photon cross sections at hadron 
colliders is sensitive (for both types of isolation cri- 
teria) on the photon fragmentation function at large 
momentum transfer, which has up to now been mea- 
sured only at LEP^^'^^. Further information on the 
photon fragmentation function at large momentum 
transfer might be gained from yet unanalyzed LEP 
data or from the study of photon-plus-jct final states 
in deep inelastic scattering at HERA'''^, where first 
data are now becoming available^^. 

4-2 Photon Pairs 

One of the most promising channels for the discovery 
of a light Higgs boson (mn < 140 GeV) at the LHC 
is based largely on the observation of the rare decay 
to two photons. To perform an accurate background 
subtraction for this observable, one requires a precise 
prediction for QCD reactions yielding di-photon fi- 
nal states. At first sight, the 0{a^) process qq 77 
yields the leading contribution. However, due to the 
large gluon luminosity at the LHC, both qg (777 
{0{al)) and gg 77 (©(a^)) subproccsscs yield 
contributions of comparable magnitude. The NLO 
corrections to the qq and qg subprocesses have been 
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Figure 6. Photon pair production for different isolation 
criteria^'' . 



ments, where one expects the main production chan- 
nel to be gluon fusion, mediated through a top quark 
loop. To a good approximation''^, one can use an ef- 
fective gluon-gluon-Higgs coupling to describe this 
process in perturbative QCD (provided the leading 
order mass dependence is factored out explicitly) . In 
this approximation, the calculation higher order cor- 
rections to inclusive Higgs production becomes very 
similar to the analogous calculation for gauge boson 
production. 

Inclusive vector boson production has been com- 
puted to NNLO^" already more than ten years ago. 
Very recently, these results have been verified for the 
first time in an independent calculation*^, carried 
out in the context of the derivation of NNLO correc- 
tions to inclusive Higgs boson production. 



known for quite some time, these are implemented in 
the flexible parton level event generator DIPHOX''^ . 
Most recently, NLO corrections were also derived for 
the gg subprocess^^. Since the lowest order contri- 
bution to this process is already mediated by a quark 
loop, this calculation contains some of the features 
appearing in jet physics only at NNLO, such as two- 
loop QCD amplitudes and unresolved limits of one- 
loop amplitudes (see Section 3.3 above). Another 
important new result are the NLO corrections to 
two-photon-plus-jet production'^^ , forming the back- 
ground to Higgs boson detection at large transverse 
momenta. 

It must be kept in mind that the di-photon cross 
sections are highly sensitive on the isolation criteria 
applied to the photons, Figure 6, with a substantial 
contribution arising from photon fragmentation at 
large momentum transfers^^. Moreover, it is experi- 
mentally difficult to distinguish photons from highly 
energetic neutral pions which decay into a closely 
collimated photon pair, mimicking a single photon 
signature. The pion background in photon pair pro- 
duction has been studied to NLO^* and implemented 
in DIPHOX, showing that in particular the 7r°7 chan- 
nel remains comparable to the 77 channel even for 
tight isolation criteria. 

5 Higgs and Gauge Boson Production 

The search for the Higgs boson is one of primary 
goals of present and future hadron collider experi- 



5.1 Higgs Boson 

The NNLO corrections to the Higgs production 
cross section were obtained first in the soft/coUinear 
approximation*^; shortly thereafter, the full coeffi- 
cient functions were obtained by expansion around 
the soft limit*^, and fully analytically*^ by extend- 
ing the IBP /LI reduction method and the differential 
equation technique (see Section 3.3) to compute dou- 
ble real emission contributions. These results were 
confirmed independently*^ using the techniques of 
the original vector boson calculation. It turned out 
that inclusion of NNLO corrections yields a sizable 
enhancement of the Higgs production cross section. 
Figure 7, and a reduction of the uncertainty due to 
renormalization and factorization scale. Recently, 
this calculation was further improved by the inclu- 
sion of effects due to soft gluon resummation*^ . Fur- 
ther NNLO results on inclusive Higgs boson produc- 
tion involve pseudoscalar Higgs production through 
gluon fusion*^, Higgs production in bottom quark 
fusion*^ and Higgs-strahlung off a vector boson**. 

Since hadron collider experiments only cover a 
limited range of the final state phase space, it is very 
desirable to have not only predictions for the inclu- 
sive Higgs production cross sections, but also for dif- 
ferential distributions in rapidity and transverse mo- 
mentum. Next-to-leading order corrections to both 
distributions*^^^"^^^ became available recently. The 
rapidity distribution. Figure 8, is only moderately 
modified, but extends out significantly beyond the 
experimental coverage. The calculation of correc- 
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Figure 7. Inclusive Higgs production at the LHC^'^. 



tions to transverse momentum distributions is reli- 
able only for transverse momenta larger or equal to 
the Higgs boson mass, while multiple soft gluon ra- 
diation plays a crucial role at lower transverse mo- 
menta. The resummation of these corrections was 
performed recently^^ to NNLL accuracy, Figure 9. 

5.2 Gauge Bosons 

The production cross sections for and bosons 
at hadron colliders are well understood both exper- 
imentally and theoretically. At present, these cross 
sections are measured to an error of about 10% from 
Tevatron Run 1, largely limited by statistics. A con- 
siderable reduction of the experimental error is an- 
ticipated from Run 11 and for the LHC. On the the- 
ory side, inclusive vector boson production has been 
computed to NNLO®°'^^ 

Given the good theoretical and experimental un- 
derstanding of and boson production, it has 
been suggested to use these for a determination of the 
LHC luminosity^^ . In practice, it turns out that it 
is not possible to measure the fully inclusive produc- 
tion cross sections, but only cross sections integrated 
over a restricted range in rapidity, for which NNLO 
corrections were also computed very recently^^. For 
the production, which is observed only through 
the Iv decay channel, it is moreover mandatory to 
compute the spatial distribution of the decay prod- 
ucts, which is however only known to NLO^^. 

A crucial ingredient to precision NNLO predic- 
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LHC, including soft gluon resummation''^. 



tions of these cross sections at LHC are parton distri- 
butions accurate to this order. The determination of 
parton distributions from a global fit to experimental 
data is described in great detail in Robert Thome's 
talk at this conference^^ . To perform a fit at NNLO, 
one needs on the one hand the coefficient functions 
for all contributing observables (deep inelastic scat- 
tering, Drell-Yan process, jet production and possi- 
bly direct photon production) to NNLO. At present, 
only Drell-Yan process and deep inelastic scattering 
are known to this order. On the other hand, also 
the partonic splitting kernels (Altarelli-Parisi split- 
ting functions) are required to NNLO accuracy. At 
the moment, this calculation is ongoing. The method 
applied in this calculation is the determination of the 
splitting functions from the forward photon-parton 
scattering amplitude at three loops, evaluated in mo- 
ments of the external partonic momentum. Interme- 
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diate results involve some fixed moments®^, as well 
as all moments for the non-singlet fermionic loops®'^. 

6 Conclusions and Outlook 

QCD affects all observables studied at present and 
future hadron colliders. Given the anticipated lumi- 
nosity of these machines, QCD reactions there will be 
precision physics, very much like electroweak physics 
was precision physics in the LEP era. The study of 
many of the standard scattering reactions will al- 
low a precise determination of the strong coupling 
constant, electroweak parameters, quark masses and 
parton distribution functions. In turn, this infor- 
mation translates in improved predictions for new 
physics signals and their backgrounds. 

Both the precision determination of standard 
model parameters (and auxiliary quantities) and the 
design of search strategies for new physics effects re- 
quire substantial input from theoretical calculations. 
For the purpose of drafting searches for physics be- 
yond the standard model, one often requires pre- 
dictions for multiparticle final states, which are at 
present only available at the leading order. Such ac- 
curacy is in general sufficient for this purpose, given 
that the available generic programs also provide in- 
terfaces to partonic showers and hadronization mod- 
els, thus predicting fully hadronic events, which can 
be further processed through detector simulations. 
Next-to-leading order calculations will be important 
to refine searches, and to identify potential new sig- 
nals, since quantitative predictions start to become 
reliable only at this order. The current frontier of 
NLO calculations are 2^3 reactions, where some of 
the most prominent observables are known. Exten- 
sion to 2 ^ 4 reactions will require new theoretical 
tools, in particular towards generic, process indepen- 
dent algorithms. Another important development in 
NLO calculations is the interface to partonic showers, 
which has recently been devised. For the precise ex- 
traction of standard model parameters from bench- 
mark reactions, NNLO calculations will be manda- 
tory, since these reactions are (already at present 
colliders) measured to a level of accuracy at which 
the theoretical error on the NLO calculation becomes 
the dominant source of uncertainty. Presently, first 
NNLO calculations were performed for 2^1 reac- 
tions, and 2^2 calculations are well under way. 
For applications at hadron colliders, NNLO calcula- 



tions become only meaningful if augmented by par- 
ton distributions accurate to NNLO, which require 
the knowledge on the three-loop splitting functions, 
which are also calculated at present. 

Many observables do moreover require the re- 
summation of large logarithms spoiling the conver- 
gence of the pcrturbativc series. Fragmentation ef- 
fects enter many observables with identified particles 
in the final state. In particular, a consistent treat- 
ment of heavy quark fragmentation effects can ac- 
count for a large part of the observed discrepancy 
in B hadron spectra, and quark-to-photon as well as 
quark-to-pion fragmentation yield important contri- 
butions to photon pair final states forming an impor- 
tant background to Higgs searches. Much valuable 
information on these fragmentation functions is con- 
tained in data from LEP, and should be extracted 
(as long as this is still a feasible task) to improve 
predictions for collider observables. 
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DISCUSSION 

Sungwon Lee (Texas A & M University): Cur- 
rently there are large discrepancies between data 
and NLO QCD predictions for direct photon 
production. Has there been any theoretical 
progress on the intrinsic (effective parton) kt is- 
sue for direct photon production? 

Thomas Gehrmann: Understanding of kt effects 
is one motivation for doing NNLO calculations. 
At NNLO, for the first time, we start to fully 
model kt effects due to hard parton emission 
in the initial state because we allow for either 
double emission of one of the incoming legs, or 
for double uncorrelated emission. So, once you 
have NNLO calculations for 2 ^ 2 scattering 
processes available, you will really have a the- 
oretical tool for computing kt effects from per- 
turbation theory. 

Ikaros Bigi (Notre Dame University): You showed 
this very instructive curve about b quark frag- 
mentation, where you showed that the Peterson 
et al. prediction is less than optimal. Do you 
have (or does someone have) a similar curve for 
charm fragmentation? 

ThomEis Gehrmann: The recent work by Cacciari 
and Nason addresses this issue. They refit 
charm fragmentation functions using a new pa- 
rameterization in moment space trying to ex- 
pose the information content of LEP data rel- 
evant to proton- antiproton colliders. Compari- 
son with old parameterizations is however not 
made. The basic message is that when you 
start fitting fragmentation functions, you have 
to be extremely careful that you are not intro- 
ducing artifacts from the choice of parameteri- 
zation which then, although giving a you a least 
fit, do not really reproduce the data because 
you are starting with too stiff an initial form. 



